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Abstract. MicroRNAs (miRNAs) are small, non-coding 
RNAs that negatively regulate gene expression and have 
emerged as potential biomarkers in radiation response to 
human cancer. Only a few miRNAs have been identified in 
radiation response to prostate cancer and the involvement of 
the radiation-associated miRNA machinery in the response of 
prostate cancer cells to radiation is not thoroughly understood. 
Therefore, the purpose of the present study was to compre- 
hensively investigate the expression levels, arm selection 
preference and isomiRs of radiation-response miRNAs in 
radiation-treated PC3 cells using a next-generation sequencing 
(NGS) approach. Our data revealed that the arm selection 
preference and 3' modification of miRNAs may be altered 
in prostate cancer after radiation exposure. In addition, the 
proportion of AA dinucleotide modifications at the end of 
the read gradually increased in a time-dependent manner 
after PC3 radiation treatment. We also identified 6 miRNAs 
whose expression increased and 16 miRNAs whose expres- 
sion decreased after exposure to 10 Gy of radiation. A pathway 
enrichment analysis revealed that the target genes of these 
radiation-induced miRNAs significantly co-modulated the 
radiation response pathway, including the mitogen-activated 
protein kinase (MAPK), Wnt, transforming growth factor-p 
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(TGF-p) and ErbB signaling pathways. Furthermore, analysis 
of The Cancer Genome Atlas (TCGA) database revealed 
that the expression of these radiation-induced miRNAs 
was frequently dysregulated in prostate cancer. Our study 
identified radiation-induced miRNA candidates which may 
contribute to radiosensitivity and can be used as biomarkers 
for radiotherapy. 

Introduction 

Prostate carcinoma is the most frequently diagnosed visceral 
cancer in men worldwide. An increasing prevalence has been 
reported in recent decades (1). Radiation therapy is one of the 
primary modalities in prostate cancer treatment. Ionizing radi- 
ation damages cells through free radicals from the radiolysis 
of water that cause DNA double-strand breaks. However, the 
efficacy of the radiotherapy may be affected by the cellular 
response to radiation. Radiotherapy is highly effective in 
treating radiosensitive tumors and enhancing the therapeutic 
efficacy can increase the overall survival rate. However, the 
presence of radioresistant tumors leads to cancer relapse and 
metastasis. Understanding the tumor-radiation-related genes 
to predict the tumor response to radiotherapy may potentially 
modulate the treatment outcome for prostate cancer patients. 

MicroRNAs (miRNAs) are a family of small, non-coding, 
single-stranded RNAs composed of ~22 nucleotides (nt) that 
negatively regulate protein expression at the post-transcrip- 
tional level (2). They function as gene regulators by binding 
to partially complementary sites of mRNAs and cause transla- 
tion inhibition or direct degradation of the target mRNA. It has 
been suggested that miRNAs are responsible for controlling 
~50% of all protein-coding genes (3). The widespread regula- 
tion of protein levels has been studied in cellular models (4). 
Previous studies have demonstrated that the expression of 
miRNAs is clearly involved in cancer development, and the 
deregulation of several miRNAs has been observed in various 
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types of cancer, including prostate cancer. Porkka et al (5) was 
the first to identify a miRNA signature specific for prostate 
cancer by systematically profiling prostate cancer cell lines. 
Numerous studies have identified many dysfunctional miRNAs 
by using a high-throughput approach, which contributed to 
prostate cancer progression, including the let-7 family, miR-1, 
-20a, -21,-34a, -106b, -125b, -205 and -521 (6-13). Although 
several studies have investigated the role of these dysfunc- 
tional miRNAs to develop prostate cancer therapy, few studies 
have determined the roles of miRNAs in radiation response 
in prostate cancer. The upregulation of miR-521 reduces 
the response to radiation damage by specifically targeting a 
DNA repair protein, the Cockayne syndrome protein A (13). 
Li et al (14) found that miR-106b was dysregulated after 
radiation treatment and suppressed radiation-induced p21 
activation, suggesting it may override radiation-induced cell 
cycle arrest and cell growth inhibition. Radiation delivered in 
daily fractions altered a greater number of miRNAs compared 
with single-dose radiation, and involved the upregulation of 
miR-34a and let-7 miRNAs (15). 

Next-generation sequencing (NGS) is a high-throughput 
screening technology, and NGS data can be applied in investi- 
gating miRNA expression, miRNA isoforms (isomiRs) and the 
arm selection preferences of miRNAs. Therefore, the purpose 
of the present study was to comprehensively investigate the 
distribution of miRNAs after radiation treatment in PC3 cells 
by using an NGS approach. Furthermore, we explored the 
function of radiation-associated miRNA by conducting an 
in silico analysis. 

Materials and methods 

Cell culture and radiation treatment. A PC3 cell line was 
obtained from the American Type Culture Collection and was 
maintained in RPMI-1640 and supplemented with 10% inac- 
tivated fetal bovine serum (FBS; Invitrogen, Carlsbad, CA, 
USA). The cells were exposed to various radiation dosages (0, 
2, 6, 10, 14 and 18 Gy) and were subsequently cultured in fresh 
medium. The total RNA was obtained at various time points 
(0, 5, 15 and 40 h after treatment) by using TRIzol (Invitrogen) 
according to the manufacturer's instructions. The concentra- 
tion, purity and amount of total RNA were determined using a 
NanoDrop 1000 spectrophotometer (NanoDrop Technologies, 
Inc., USA). 

Collection and preprocessing of sequence reads. PC3 cells 
were exposed to 10 Gy of radiation. After radiation treatment, 
the cells were lysed at various time points (0, 5, 15 and 40 h) 
for RNA extraction. The RNA samples were prepared using 
an Illumina small RNA preparation kit, and were subsequently 
sequenced using the Illumina HiSeq platform. The generated 
sequence reads were first subjected to quality control to remove 
low-quality reads. The sequence reads were then subjected to 
3' adaptor trimming to generate clean reads, as previously 
described (15,16). To attain a high confidence level, only the 
clean reads with a read count >2 and with a length ranging 
from 15 to 27 nt were included in further analyses. 

Mapping clean reads to pre-miRNAs. To investigate miRNA 
expression profiles in different libraries, we mapped the qualified 



clean reads back to human pre-miRNAs (miRBase 19). To 
eliminate ambiguous multiple hits during the mapping proce- 
dure, no mismatch was allowed. Previous studies reported that, 
when mapped back to pre-miRNAs, sequence reads usually 
carried mismatches preferentially located at their terminal 3' 
ends (17-20). This mismatch was named the 3' end modification. 
To determine whether the 3' end modification patterns differed 
among libraries, as described in our previous studies (21), we 
trimmed and collected the terminal 3' end mismatches one by 
one. In addition, the remaining perfect match reads had to be 
at least 18 nt in length. As a result, we kept reads with no less 
than 18-nt perfect alignment and 3' end modification patterns. 

Classifying non-miRNA reads into different data sets. The 
sequence reads that may not be mapped back to pre-miRNAs 
were classified into classes by mapping to acquire different 
data sets with Bowtie (22) and allowing a single nucleotide 
variation. The sequences of mRNAs and other ncRNAs were 
derived from the NCBI RefSeq 47 (23). The tRNA sequences 
were downloaded from the Genomic tRNA database (24) and 
the rRNA sequences were downloaded from the SILVA data- 
base (25). The snoRNA, scaRNA and snRNA sequences were 
all downloaded from NONCODE (26). The sequence reads not 
belonging to any of the described RNA classes were uploaded 
to the RepeatMasker to identify repeat elements, which were 
classified as unknown. 

miRNA expression level according to The Cancer Genome 
Atlas (TCGA) data. TCGA project collects both cancer and 
corresponding normal tissues from hundreds of prostate 
cancer patients. We downloaded all level-3 miRNA expres- 
sion data of prostate adenocarcinoma from the TCGA Data 
Portal (https://tcga-data.nci.nih.gov/tcga/dataAccessMatrix. 
htm). These level-3 data included calculated expressions for 
each miRNA derived from the Illumina HiSeq sequencing 
results. A total of 198 tumor samples and 50 normal samples 
were found at the time the data were downloaded. We kept 
only the expression data of 50 participants who had both 
miRNA expression levels from both tumor and normal 
tissues. Normalized quantification expression levels for these 
50 participants were further examined for each investigated 
miRNA. 

Pathway enrichment analysis. We attempted to determine 
the functions of the miRNA target genes by investigating 
the pathways with which the miRNA target genes were 
involved. Therefore, we first downloaded the target genes of 
differentially expressed miRNAs from TargetScan 6.0, and 
then mapped the target genes onto the Kyoto Encyclopedia 
of Genes and Genomes (KEGG) pathways based on the 
Enzyme Commission (EC) numbers by using the R package 
SubPathwayMiner v.3.1 (27). Subsequently, the hypergeo- 
metric test was performed to identify significantly enriched 
pathways and calculate the false positive discovery rate in the 
FDR-corrected q-value. 

Results 

miRNA profiling of radiation-treated prostate cancer cells. 
To characterize the mechanism involved in the radiation 
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Table I. Summary of sequence reads and the detected miRNAs. 
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Figure 1. Radiation treatment of human prostate cancer cells, PC3. (A) PC3 
cells were treated with various radiation doses (0, 2, 4, 6, 8, 10, 12, 14, 
16 and 18 Gy) and were subsequently subjected to fresh culture medium. 
After culturing for an additional 4 days, the morphology was observed using 
light microscopy (x40 magnification). (B) The expression pattern of COX-2 
and p21 in radiation-treated PC3 cells was examined using a real-time PCR 
method. S26 was used as an internal control. (C) The distribution of small 
RNA reads in 11 categories was classified. 



response of prostate cancer, we used NGS to comprehensively 
analyze the distribution of miRNAs after radiation treat- 
ment in PC3 cells. As indicated in Fig. 1A, PC3 cells were 
exposed to various dosages of radiation (0, 2, 4, 6, 8, 10, 
12, 14, 16 and 18 Gy) and then subjected to a fresh culture 
medium for an additional 4 days. We found that the growth of 
the PC3 cells obviously decreased when exposed to 10 Gy of 
radiation. Therefore, we collected cell RNA at various times 
(0, 5, 15 and 40 h) following the 10-Gy radiation treatment. 
We confirmed the expression levels of Cox-2 and p21, which 
may be induced by radiation at 24 h according to previous 
studies (16,28). The expression levels of Cox-2 and p21 may 
be upregulated by radiation treatment in PC3 cells (Fig. IB). 
We then performed the comprehensive miRNA profile at 
various time-points in radiation-treated PC3 cells by using the 
Illumina HiSeq platform. 
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712 
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10,589,934 


80.86 


739 
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P_C, P_5, P_15 and P_40 are prostate cancer cell lines with different 
radiation treatment. By mapping the clean sequence reads back to 
pre- miRNAs, we can quantify how many pre- miRNAs and mature 
miRNAs were detected. 



Table II. miRNAs with altered expression in response to radia- 
tion in PC3 cells using next-generation sequencing. 
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hsa-miR-9-5p 


1 


2.82 


0.90 


2.63 




hsa-miR-22-3p 


1 


1.60 


2.85 


2.73 




hsa-miR-25-3p 


1 


1.54 


1.39 


2.33 


Upregulation d 


hsa-miR-30a-5p 


1 


2.81 


2.43 


2.33 


Upregulation 0 


hsa-miR-550a-3p 


1 


1.88 


1.46 


2.09 


Upregulation d 


hsa-miR-548h-5p 


1 


0.50 


0.44 


2.56 




Downregulation b 












hsa-let-7c 


1 


0.93 


0.78 


0.45 




hsa-let-7d-5p 


1 


0.30 


0.11 


0.40 


Upregulation d 


hsa-let-7e-5p 


1 


0.58 


0.55 


0.40 




hsa-miR-15a-5p 


1 


0.78 


0.67 


0.45 


Upregulation' 1 


hsa-miR-17-3p 


1 


0.52 


0.49 


0.47 


Upregulation' 1 


hsa-miR-30d-3p 


1 


0.92 


0.75 


0.41 


Upregulation' 1 


hsa-miR-92a-5p 


1 


0.65 


0.52 


0.50 


Upregulation' 1 


hsa-miR-125a-3p 


1 


0.42 


0.42 


0.32 




hsa-miR-197-3p 


1 


0.77 


0.79 


0.44 


Upregulation' 1 


hsa-miR-221-5p 


1 


0.59 


0.31 


0.44 


Downregulation d 


hsa-miR-320b 


1 


0.65 


0.41 


0.32 


Upregulation* 1 


hsa-miR-342-5p 


1 


0.59 


0.63 


0.47 


Upregulation* 1 


hsa-miR-361-3p 


1 


0.45 


0.53 


0.40 


Upregulation* 1 


hsa-miR-374a-5p 


1 


1.04 


0.93 


0.47 




hsa-miR-501-3p 


1 


0.77 


0.48 


0.45 


Upregulation* 1 


hsa-miR-671-3p 


1 


0.77 


0.62 


0.41 


Upregulation* 1 



"Expression levels of miRNA were inducted >2-fold change after PC3 
radiation treatment with 10 Gy for 40 h. b Expression levels of miRNA 
were repressed >2-fold change after PC3 radiation treatment with 
10 Gy for 40 h. c d The difference was indicated to be significant with 
p- value less than 0.01 or 0.001 . 



Analysis of miRNA sequence reads. Once the samples 
were sequenced, we collected >9 million clean reads in all 
libraries (Table I). In addition to miRNA, we also determined 
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Figure 2. The distribution of isomiRs and -5p/-3p arm selection in PC3 radiation treatment. (A) Mapping results of hsa-miR-2110. As shown in 5p:8-29, hsa- 
miR-21 10 encodes mature miRNA at only its 5 p arm, and miRNA spans from nucleotide 8 to nucleotide 29 of the hairpin. The integer values on the left denote 
the read count of each isomiR. The comma-separated values denote the position shift in the isomiR relative to the miRBase annotated positions (8 to 29). The 
nucleotides in lowercase type denote the sequence fragments originating from the 3' modification event. (B) The proportion of 3' modifications at end of reads 
at different time-points after PC3 radiation treatment as observed from NGS data. (C) Fold-change of the -5p/-3p arm of miRNA at different time points after 
PC3 radiation treatment as observed from NGS data. isomiRs, miRNA isoforms; miRNA, microRNA; NGS, next-generation sequencing. 



which molecules were the remaining non-miRNA reads. By 
mapping the non-miRNA reads back to a different data set, 
we classified the reads into 11 categories. Fig. 1C demonstrates 
that miRNA accounted for 80% of all clean reads in the pros- 
tate cell libraries. Other categories accounted for relatively 
low proportions, which indicated the high performance of the 
sample preparation protocol. In addition, the proportions of 
the categories were considerably similar among libraries, indi- 
cating that radiation treatment did not alter the composition of 
RNA samples in the prostate cell libraries. 

After mapping the clean reads to the genome, most of 
the miRNA reads tend to exist as isomiR. As demonstrated 
in Fig. 2A, hsa-miR-2110-5p had 5 isomers, whereas the 
opposite-arm miRNA-3p had 15 isomiRs, which demonstrated 
that abundant miRNAs tend to have more isomiRs. Our data 
revealed that the isomiR quantity was highly correlated with 



miRNA abundance (Pearson's correlation coefficient, 0.91). 
In addition, we observed that the modified nucleotides were 
preferentially located at the 3' end of the sequence read 
(presented in lower case in Fig. 2A). The data indicated that 
one A nucleotide or one U nucleotide was frequently added at 
the end of the read. Notably, we found that the proportion of 
AA dinucleotides modified at the end of the read was gradu- 
ally increased in a time-dependent manner after the PC3 cells 
were treated with radiation, which indicated that the 3' end 
modification may be altered by radiation treatment in PC3 
cells. Our previous studies indicated that the use of miR-5p 
and -3p may be altered in human cancer (29-31). In the present 
study, our data indicated that arm selection preference was 
consistent across nearly all libraries. Only a few cases were 
observed in which the use of -5p and -3p arm selection had 
different preferences at various time-points after radiation 
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Figure 3. (A) Expression levels of radiation-induced miRNAs in prostate cancer. (B) Expression levels of radiation-suppressed miRNAs in prostate cancer. 
Expression levels of miRNAs between tumor and corresponding normal tissues from 50 prostate cancer patients were analyzed using TCGA dataset. The 
expression levels of miRNAs were presented in transcript per million (TPM). The expression level between tumor and normal cells was evaluated by conducting 
paired t-tests (P<0.05 was considered significant; NS, non-significant. *P<0.05, ™P<0.01, *™P<0.001). TCGA, The Cancer Genome Atlas. 



treatment (Fig. 2C). Further research is required to support 
these findings. 

Radiation-response miRNAs in prostate cancer. By summarizing 
the read count of all the isomiRs that belonged to the same mature 
miRNAs, we quantified the miRNA expression abundances, and 
presented the result in transcript per million (TPM). Twenty-two 
miRNAs were selected and are presented in Table II, demon- 
strating that their expression levels were altered >2-fold after 
being subjected to radiation exposure (expression of 6 miRNAs 
increased, and expression of 16 miRNAs decreased). To explore 
the putative role contributing to prostate cancer progression, we 
examined the effect of expression levels of radiation-associated 
miRNAs on prostate cancer from the available TCGA dataset by 
using an in silico analysis. We downloaded 100 miRNA expres- 
sion profiles from 50 prostate cancer patients, including 50 cancer 
lesion and 50 corresponding normal tissues. As demonstrated 
in Fig. 3A, the expression levels of radiation-induced miRNAs, 
miR-25, miR-30a and miR-550a, were significantly upregulated 
in the prostate cancer cells compared with the corresponding 
normal tissue cells. Twelve radiation-suppressed miRNAs were 
identified, i.e. let-7d, miR-15a, miR-17, miR-30d, miR-92a, 
miR-197, miR-221, miR-320b, miR-342, miR-361, miR-501 
and miR-671, and a significantly different expression between 
prostate cancer and the corresponding adjacent part was found, 
including 11 upregulated and 1 downregulated (Fig. 3B). Overall, 
the data indicated that most of the radiation-response miRNAs 
were identified as dysregulated in prostate cancer according to 
an in silico analysis (15/22; 1 downregulated, 14 upregulated and 
the rest demonstrated no change in expression in prostate cancer). 

Pathway enrichment analysis of miRNAs. miRNAs can func- 
tion as either oncogenes or tumor suppressors depending on 
their target genes. Therefore, identifying a target can facilitate 



elucidating the role of miRNAs in prostate cancer treatment 
radiation (32,33). Typically, one miRNA tends to have hundreds 
of target genes and a group of miRNAs co-modulated as a 
biological function involved in the regulation of a signaling 
pathway. Therefore, we further explored the biological function 
of radiation-response miRNAs by conducting a pathway- 
enrichment analysis. The putative target genes of miRNAs 
were obtained from TargetScan 6.0; subsequently, these target 
genes of the individual miRNAs were mapped onto KEGG 
pathways. Our data indicated that the target genes of radiation- 
response miRNAs were frequently significantly enriched in 
several cancer- or radiation-related pathways, including the 
mitogen- activated protein kinase (MAPK), ErbB, p53, Wnt, 
transforming growth factor-p (TGF-p) and mTOR signaling 
pathways with an FDR <0.05 (Table III). We also subjected 
the target genes of the 2-gene set, upregulated miRNA and 
downregulated miRNAs, to pathway-enrichment analysis. 
Similar results were observed; their targets were significantly 
enriched in the prostate cancer pathway (Fig. 4) and radiation- 
related pathways, including the MAPK, ErbB, Wnt and TGF-P 
signaling pathways (Tables IV and V). 

Discussion 

Our previous studies indicated that the distributions of 3' end 
modifications and the arm selection preference of miRNAs 
were different between normal and tumor tissues (29-31). 
The -5p and -3p of miRNA play a distant role by suppressing 
the different target genes. It was previously reported that, in 
contrast to the oncogenic effect of miR-17 (-5p), miR-17*(-3P) 
plays a tumor suppressive role in prostate cancer (9,34,35). The 
miR-28-5p and miR-28-3p also play opposite roles in colon 
cancer cell proliferation and migration (36). In the present 
study, our data showed that the -5p and -3p of particular 
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Table III. The enriched pathways of radiation-induced miRNA target genes. 


microRNA 


Cancer-relative pathway (FDR<0.05) 


Upregulation 




hsa-miR-9-5p 


Focal adhesion, pathways in cancer. ErbB signaling pathway. MAPK signaling pathway. 




prostate cancer 


hsa-miR-22-3p 


Chronic myeloid leukemia. MAPK signaling pathway. ErbB signaling pathway, pathways 




in cancer glioma nrostate cancer nhosnhatidvlinositol si^nalin? system colorectal cancer 


hsa-miR-25-3p 


N.D 


Uafl llliA Jud JU 


n n 


IlSa-IIlUV-JJUa-jp 




nsa-miKo^fonop 


IN .YJ 


Downregulation 




hsa-let-7c 


MAPK signaling pathway, pathways in cancer. p53 signaling pathway, melanoma. 




chronic myeloid leukemia, glioma, pancreatic cancer, focal adhesion, small cell lung cancer, 




bladder cancer, prostate cancer 


hsa-let-7d-5p 


MAPK signaling pathway, pathways in cancer. p53 signaling pathway, melanoma. 




chronic myeloid leukemia, glioma, pancreatic cancer, focal adhesion, small cell lung cancer, 




bladder cancer, prostate cancer 


hsa-let-7e-5p 


MAPK signaling pathway, pathways in cancer. p53 signaling pathway, melanoma. 




chronic myeloid leukemia, glioma, pancreatic cancer, focal adhesion, small cell lung cancer, 




bladder cancer, prostate cancer 


hsa-miR-15a-5p 


Pathways in cancer, regulation of actin cvtoskeleton. renal cell carcinoma. MAPK signaling pathway. 




focal adhesion, melanoma, prostate cancer. Wnt signaling pathway. p53 signaling pathway. 




mTOR signaling pathway, non-small cell lune cancer, pancreatic cancer, cell cycle 


hsa-miR-17-3p 


MAPK signaling pathway, pathways in cancer, chronic myeloid leukemia, pancreatic cancer. 




melanoma, bladder cancer. TGF-p signaling pathway, prostate cancer. mTOR signaling pathway. 




non-small cell lung cancer, renal cell carcinoma, cell cycle, p53 signaling pathway 


hsa-miR-30d-3p 


N.D 


hsa-miR-92a-5p 


N.D 


hsa-miR-125a-3p 


MAPK signaling pathway, adherens junction, pancreatic cancer. TGF-p signaling pathway 


hsa-miR-197-3p 


N.D 


hsa-miR-221-5p 


Wnt signaling pathway. ErbB signaling pathway 


hsa-miR-320b 


Chronic myeloid leukemia, non-small cell lung cancer, glioma, pathways in cancer, focal adhesion, 




pancreatic cancer, melanoma. ErbB signaling pathway, colorectal cancer. TGF-p signaling pathway. 




prostate cancer. MAPK signaling pathway. mTOR signaling pathway 


hsa-miR-342-5p 


N.D 


hsa-miR-36 1 -3p 


Pathways in cancer. mTOR signaling pathway, melanoeenesis. renal cell carcinoma. 




nucleotide excision repair 


hsa-miR-374a-5p 


Pathways in cancer, prostate cancer. TGF-p signaling pathway, endometrial cancer. 




non-small cell lune cancer, basal cell carcinoma. MAPK signaling pathway 


hsa-miR-501-3p 


N.D 


hsa-miR-671-3p 


N.D 



miRNAs were differently regulated by radiation (shown 
in Fig. 2C). Several studies have demonstrated that miRNAs 
contain various ends, which were caused by either RNA 
editing or non-template nucleotide additions (17,18,20). These 
miRNA isoforms (isomiRs) contribute to increased miRNA 
stability or strengthened miRNA-target gene interaction and 
are differentially expressed in different cellular conditions, 
including cancer (16,37,38). Our data revealed that the propor- 
tion of AA dinucleotide modifications at the end of the read 
gradually increased in a time-dependent manner after the PC3 
cells were treated with radiation, suggesting that radiation 



may influence the particular miRNA stability or efficiency of 
silencing targets by regulating the 3' end modifications, which 
warrants further research. 

miRNAs are known to function as gene silencers and 
are involved in modulating biological functions, including 
cell growth, apoptosis, the cell cycle and the metastasis of 
cancer (39). Comprehensive miRNA profiling of prostate 
cancer has indicated that several miRNAs are differentially 
expressed between prostate cancer and the adjacent normal, 
which contributes to prostate cancer progression (40-42). 
In the present study, we analyzed miRNA expression from 
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Figure 4. The enriched pathway of the target gene union of radiation-response miRNAs. (A) The target gene union of radiation-upregulated miRNAs enriched 
in the prostate cancer pathway (FDR=0.001). (B) The target gene union of radiation-downregulated miRNAs enriched in the prostate cancer pathway 
(FDR=3 .6E-9). The target genes are labeled in red. 



1 074 LEUNG et at RADIATION-INDUCED microRNA EXPRESSION IN PROSTATE CANCER 



Table IV. The enriched pathways of radiation-upregulated miRNA target genes. 



Pathway Id 


Pathway name 


FDR 


Path:04360 


Axon guidance 


H OOT7 AO 

7.88E-08 


rath:04 1 11 


Neurotrophin signaling pathway 


1 1 1 17 r\H 


n *.i a a a 1 a 

Path:04010 


MAPK signaling pathwav 


1.11E-07 


n^4-U .A CO A A 

rath:05200 


Pathways in cancer 


3.25b-0/ 


Path:04012 


ErbB signaling pathwav 


3.14E-06 


Path:04l20 


Ubiquitin mediated proteolysis 


a oo l — 1 i\ s 

8.82E-06 


"Tl 4.1 A A 1 A A 

Path:04144 


Endocytosis 


1 .24E-05 


Path:04520 


Adherens junction 


5.86E-05 


T* ' 1 A C A 1 r\ 

Path:054l2 


Arrhythmogenic right ventricular cardiomyopathy (ARVC) 


7.26E-05 


Fatn:U4MU 


Focal adhesion 


A AAA 1 AC 

U.UUU1U5 


Path:04916 


Melanogenesis 


A A A AO A C 

0 .000245 


Fatn:U5ZZU 


Chronic myeloid leukemia 


A AAA") A C 


ratn.U4olU 


Regulation of actin cytoskeleton 


U .UUUZo 1 


ratn.lozl4 


Glioma 


U. UUUZO 1 


Path:05414 


Dilated cardiomyopathy 


A A A A/^ S~ O 

0 .000268 


Path:04910 


Insulin signaling pathway 


A A A A f~\ C} s~\ 

0.000282 


Fath:04/20 


Long-term potentiation 


U .UUUZ94 


Path:04070 


Phosphatidvlinositol signaling svstem 


a AAA^A/I 

0 .000294 


T"-l .1 Af I 1 A 

Path:05410 


Hypertrophic cardiomyopathy (HCM) 


a AAAOAO 

0.000802 


Fath:05100 


Bacterial invasion of epithelial cells 


U.UU1U15 


n^4-u .aco 1 c 

Fath:05215 


Prostate cancer 


A A A 1 AA/- 

u.uuiuyo 


Fath:04310 


Wnt signaling pathwav 


r\ r\r\ 1 1 co 

0.001153 


n~*-u .r\ Am a 

Fath:04914 


Progesterone-mediated oocyte maturation 


A AAO C /I O 

0 .002548 


n 4.1 Af 01 1 

Path:05211 


Renal cell carcinoma 


a A AO 1 O C 

0.003135 


Path:04920 


Adrpocytokine signaling pathway 


A A A O A A/"\ 

0.003902 


r-i . I r\ A 1 c a 

Path:04350 


1 Gt -p signaling pathwav 


A AAOAAO 

0.003902 


"Tl 4.1 A A £L £L £L 

Path:04666 


Fc y R-mediated phagocytosis 


A A A A 1 1 A 

0.004314 


Tlr.tU.A/1 1/11 

ratn:U4141 


Protein processing in endoplasmic reticulum 


0.UU5 /38 


ratn:U5zlU 


Colorectal cancer 


U .UU5 /4 / 


T* ' 1 A <1 1 O A 

Path:04l30 


SNARE interactions in vesicular transport 


A A A^ -1 

0.007624 


T"-i ,i Ar^ i /■ 

Path:05216 


Thyroid cancer 


A A A o in 

0.008679 


Path:00562 


Inositol phosphate metabolism 


0.008679 


Path:00532 


Glycosaminoglycan biosynthesis-chondroitin sulfate 


0.008679 


n„ t L .ACA 1 /I 

Fatn:U5U14 


Amyotrophic lateral sclerosis (ALS) 


A AAO/C7A 

u.uuoo/y 


Path:05212 


Pancreatic cancer 


0.008679 


T-i .1 A/IA^A 

Path:04020 


Calcium signaling pathway 


A AAO£ in 

0.008679 


Path:05218 


Melanoma 


A A 1 A 1 

0.010147 


t"-i ,i a /i Aon 

Path:04930 


T 1 TT J " 1_ * 11"*. 

Type II diabetes melhtus 


A A 1 OTT7 

0.018227 


Path:04962 


Vasopressin-regulated water reabsorption 


0. 018227 


Path:04530 


Tight junction 


A A 1 AO OA 

0.019289 


n^«-u .a /i c 1 o 

Fath:04512 


ECM-receptor interaction 


a a 1 noon 
0.019289 


Path:049l2 


GnRH signaling pathway 


A A 1 AOOA 

0. 019289 


Path:05223 


Non-small cell lung cancer 


A AOOO A A 

0.023349 


Fath:05222 


Small cell lung cancer 


A AOO /IOC 

0.032425 


Path:05213 


Endometrial cancer 


A f\1 C A A 

0.035464 


Path:04662 


B cell receptor signaling pathway 


0.035464 


Path:05131 


Shigellosis 


0.035464 


Path:05221 


Acute myeloid leukemia 


0.03804 


Path:04540 


Gap junction 


0.041068 


Path:00250 


Alanine, aspartate and glutamate metabolism 


0.045237 


Path:04114 


Oocyte meiosis 


0.049286 
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Pathway Id 


Pathway name 


FDR 


Path:04010 


MAPK signaling pathwav 


0 


Path:04360 


Axon guidance 


0 


Path:05200 


Pathways in cancer 


0 


Path:04722 


Neurotrophin signaling pathway 


1.18E-12 


Path:04310 


Wnt signaling pathwav 


1.50E-10 


Path:04510 


Focal adhesion 


4.27E-10 


Path:04144 


Endocytosis 


1.40E-09 


Path:04810 


Regulation of actin cytoskeleton 


2.34E-09 


Path:05215 


Prostate cancer 


3.58E-09 


Path:05211 


Renal cell carcinoma 


3.58E-09 


Path:04720 


Long-term potentiation 


1.14E-08 


Path:05220 


Chronic myeloid leukemia 


341E-08 


Path:04120 


Ubiquitin mediated proteolysis 


740E-08 


Path:04020 


Calcium signaling pathwav 


1.10E-07 


Path:05212 


Pancreatic cancer 


1.57E-07 


Path:05214 


Glioma 


1.59E-07 


Path:05218 


Melanoma 


2.36E-07 


Path:05223 


Non-small cell lung cancer 


3.64E-07 


Path:04916 


Melanogenesis 


4.67E-07 


Path:04520 


Adherens junction 


4.67E-07 


Path:04910 


Insulin signaling pathway 


6.79E-07 


Path:04350 


TGF-P signaling pathwav 


7.72E-07 


Path:05210 


Colorectal cancer 


2.50E-06 


Path:04012 


ErbB signaling pathwav 


3.89E-06 


Path:05222 


Small cell lung cancer 


3.20E-05 


Path:04730 


Long-term depression 


5.13E-05 


Path:05221 


Acute myeloid leukemia 


7.50E-05 


Path:04150 


mTOR signaling pathwav 


8.68E-05 


Path:05213 


Endometrial cancer 


8.68E-05 


Path:05217 


Basal cell carcinoma 


9.44E-05 


Path:04710 


Circadian rhythm-mammal 


9.64E-05 


Path:04070 


Phosphatidylinositol signaling system 


9.64E-05 


Path:04540 


Gap junction 


9.83E-05 


Path:04115 


p53 signaling pathwav 


0.000113 


Path:04141 


Protein processing in endoplasmic reticulum 


0.000137 


Path:04114 


Oocyte meiosis 


0.000196 


Path:05014 


Amyotrophic lateral sclerosis (ALS) 


0.000353 


Path:04970 


Salivary secretion 


0.000414 


Path:04660 


T cell receptor signaling pathway 


0.000449 


Path:04666 


Fc y R-mediated phagocytosis 


0.000526 


Path:04512 


ECM-receptor interaction 


0.000526 


Path:05142 


Chagas disease 


0.000563 


Path:04062 


Chemokine signaling pathway 


0.000812 


Path:04210 


Apoptosis 


0.000844 


Path:04914 


Progesterone-mediated oocyte maturation 


0.001465 


Path:04110 


Cell cvcle 


0.001527 


Path:04662 


B cell receptor signaling pathway 


0.00169 


Path:04530 


Tight junction 


0.001692 


Path:04930 


Type II diabetes mellitus 


0.002128 


Path:05412 


Arrhythmogenic right ventricular cardiomyopathy (ARVC) 


0.002128 


Path:04920 


Adipocytokine signaling pathway 


0.002184 


Path:04664 


Fc 8 RI signaling pathway 


0.002621 
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Table V. Continued. 


Pathway Id 


Pathway name 


FDR 


Path:04960 


Aldosterone-regulated sodium reabsorption 


0.002802 


Path:04912 


GnRH signaling pathway 


0.002955 


T"-\ , 1 /-\ A ^\ /"\ f\ 

Path:04320 


Dorso-ventral axis formation 


0.00373 


Path:05219 


Bladder cancer 


0.00373 


Path:04340 


Hedgehog signaling pathway 


0.005218 


Path:04971 


Gastric acid secretion 


0.005633 


Path:04l30 


SNARE interactions in vesicular transport 


0.006926 


Path:00532 


Glycosaminoglycan biosynthesis-chondroitin sulfate 


0.006926 


Path:05l3l 


Shigellosis 


r\ aao 1 c>c 

0.008 185 


Path:05l60 


Hepatitis C 


f\ Ann 1 c> C 

0.008185 


Path:04630 


Jak-STAT signaling pathway 


0.008368 


Path:054l0 


Hypertrophic cardiomyopathy (HCM) 


0.009939 


Path:054l4 


Dilated cardiomyopathy 


0.010733 


Path:00534 


Glycosaminoglycan biosynthesis-heparan sulfate 


0.011359 


Path:04670 


Leukocyte transendothelial migration 


0.011852 


Path:04370 


VEGF signaling pathway 


0.013265 


Path:00562 


Inositol phosphate metabolism 


0. 013625 


Path -04270 


Va^nilar smooth mn^flp pnntrnftinn 


0 014557 


Path:00512 


O-Glycan biosynthesis 


0.016472 


Path:04330 


Notch signaling pathway 


0.026769 


Path:04142 


Lysosome 


0.038882 


Path:00533 


Glycosaminoglycan biosynthesis4ieratan sulfate 


0.038882 


Path:05145 


Toxoplasmosis 


0.048785 



TCGA database and found that the expression levels of 
radiation-induced miRNAs were frequently dysregulated in 
prostate cancer (Fig. 3). Our results are consistent with those 
of previous studies and demonstrated that miR-25, miR-17, 
miR-30d and miR-92a are overexpressed, and miR-221 is 
downregulated in prostate cancer (9,42-44). However, the 
expression levels of let-7d and miR-15a decreased according 
to TCGA, which contradicted the results of previous 
studies (45-47). These dysfunctional miRNAs have potential 
to be used as biomarkers for prostate cancer prognosis or diag- 
nosis. Therefore, understanding the function of miRNAs may 
provide practical benefits for clinical applications. Predicting 
the outcome of cancer treatment is the most promising 
application of miRNAs. Gonzales et al found miR-141 to be 
consistent with changes in other conventional biomarkers and 
to the clinical outcomes, suggesting that miR-141 can be used 
as a marker for monitoring therapeutic response in prostate 
cancer patients (48). The prognostic value of miRNA expres- 
sion profiling in prostate cancer has also been demonstrated 
by Hulf et al. They demonstrated that DNA methylation and 
histone H3K9-deacetylation of the miR-205 locus is associ- 
ated with miRNA silencing and deregulation of MED1, 
which is predictive of a poor prognosis in localized prostate 
cancer (49). 

Ionizing radiation is one of the 3 primary modalities 
used in cancer therapy. Radiation induces considerable DNA 
damages, which, if not repaired, cause cancer cells to prog- 
ress to apoptosis and cell cycle arrest. Some cancer cells are 
resistant to radiation treatment due to activation of complex 



signaling pathways that counteract these damages, including 
ErbB, nuclear factor kB (NFkB), MAPK, PI3K/AKT and trans- 
forming growth factor-p (TGF-p) signaling pathways (50-52). 
Several radiation-related miRNAs have been identified that 
contribute to the radiosensitivity of cancer cells by modu- 
lating the radiation-response signaling pathway (51,52). Since 
miRNAs are generally slightly repressed by their target 
genes, the alteration of an individual miRNA is insufficient 
for accomplishing a biological function. Previous studies 
have introduced the concept of miRNA regulatory modules 
(MRMs), which potentially serve as a model for understanding 
the detailed influences of miRNAs in cellular biological 
functions (53-55). Therefore, in the present study, we were 
particularly interested in the consequences of changes in a 
group of radiation-induced miRNAs in prostate cancer. Our 
data indicated that targets of the co-expressed miRNAs were 
enriched in a radiation-related signaling pathway, suggesting 
that they co-modulated an abundance of target genes in the 
same pathway (Table III). 

miRNAs regulate various factors in radiation-related 
biological pathways and may affect the radiosensitivity of 
tumor cells (51). Radiation-response miRNAs have been identi- 
fied in prostate cancer by using a microarray approach (13-15). 
By comparing these data, we identified known and unknown 
radiation-response miRNAs in prostate cancer by using an 
NGS approach. Li et al reported that the expression levels of 
miR-9, miR-22 and miR-30a decreased in radiation-treated 
PC3 cells (14). Radiation reduced the expression level of an 
miR-17-92a cluster and the let-7 family in prostate cancer (15). 
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In the present study, we also identified radiation-response 
miRNAs that had been reported in other types of cancer but 
not in prostate cancer, such as miR-25, miR-15a, miR-30d, 
miR-125a, miR-221 and miR-342 (21,56-63). In addition, we 
identified a group of radiation-response miRNAs that have not 
been reported in any type of cancer. The pathway-enrichment 
analysis revealed that their targets are frequently enriched in 
the radiation-response signaling pathway. 

In summary, in the present study, we thoroughly investi- 
gated radiation-response miRNAs, which may be involved in 
the radiosensitivity of prostate cancer, by modulating radiation- 
related signaling pathways using an NGS approach. These 
miRNA candidates may be effective targets for improving the 
efficacy of radiation treatment in future prostate cancer therapy. 
In addition, we observed that 3' end modifications and the -5p/-3p 
arm selection of miRNAs were altered in prostate cancer after 
radiation treatment. These finding require further research. 
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